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Concerted effort is currently ongoing to open up the Epoch of Reionization (z ~ 15-6) for stud¬ 
ies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-a 
emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct 
detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a 
direct detection is expected in the coming years, with ongoing surveys, and could open up the 
entire universe from z ~6-200 for astrophysical and cosmological studies, opening not only the 
Epoch of Reionization, but also its preceding Cosmic Dawn (z ~30-15) and possibly even the 
later phases of the Dark Ages (z ~200-30). All currently ongoing experiments attempt statistical 
detections of the 21-cm signal during the Epoch of Reionization, with limited signal-to-noise. 
Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher 
redshifts will remain out of reach. The Square Kilometre Array (SKA) will revolutionize the 
field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over 
most of the redshift range z ^6-28 with SKA 1-LOW, and possibly even higher redshifts with the 
SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy 
with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm 
emission, the different phases the universe is thought to go through, and the observables that the 
SKA can probe, referring where needed to detailed chapters in this volume. This is done within 
the framework of the current SKA1 baseline design and a nominal CD/EoR straw-man survey, 
consisting of a shallow, medium-deep and deep survey, the latter probing down to ~1 mK bright¬ 
ness temperature on arc-minute scales at the end of reionization. Possible minor modifications to 
the design of SKA1 and the upgrade to SKA2 are discussed, in addition to science that could be 
done already during roll-out when SKA1 still has limited capabilities and/or core collecting area. 
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1. Context and Layout 

In this review chapter the impact of the Square Kilometre Array (hereafter SKA1 and SKA2 for the 
first and second phases of construction) in the field of high-redshift 21-cm observations of neutral 
hydrogen during the EoR and CD is outlined, building partly on the white paper of Mellema et al. 
(2013) and referring to more than a dozen related science chapters in this volume that discuss 
particular aspects in much greater detail. In Section |] an overview of processes occurring during 
the Dark Ages, Cosmic Dawn and Epoch of Reionization (hereafter DA/CD/EoR) eras are given. In 
Section ^]the main physical processes that take place during these eras are discussed, with reference 
to the relevant accompanying chapters in this volume. In Section j| observables of the redshifted 
21-cm line are discussed, whereas in Section || we discuss the relevant parameters for any 21-cm 
survey design. In Section we discuss a three-tiered survey with the SKA of the redshifted 21-cm 
emission from the CD/EoR. In Section^ we discuss what can be done with a half (50% collecting 
area) and a full SKA1 and with SKA2 (nominally 4 times SKA1), during the rollout and build¬ 
out of SKA1 to SKA2. In Section^ we draw some general conclusions and provide suggestions. 
Throughout this chapter when referring to SKA1/2, we mean SKA1/2-LOW. 

2. A Short History of Neutral Hydrogen from Recombination to Reionization 

This section shortly summarizes the main events that occurred in the Universe between recombi¬ 
nation and reionization (see Fig.[j]), focussing on those aspects that are relevant to interferometric 
measurements of the 21-cm hypertine transition line of neutral hydrogen with the SKA. We refer to 
Barkana & Loeb (2001); Furlanetto & Briggs (2004); Lewis & Challinor (2007); Pritchard & Loeb 
(2008); Morales & Wyithe (2010); Loeb & Furlanetto (2013); Fialkov et al. (2014a) for many more 
details. 

2.1 Emission and absorption at 21-cm from z~1100 to z~6 

Just after recombination at z ~ 1100, the spin temperature of neutral hydrogen coupled very effec¬ 
tively to the neutral gas temperature, which itself coupled to the Cosmic Microwave Background 
(CMB hereafter) photons via Compton heating because of a small trace density of electrons, mak¬ 
ing it impossible to observe 21-cm radiation either in emission or absorption. This period can truly 
be regarded as an "Age of Ignorance" in cosmology with no known directly observable tracer of 
either dark matter or baryons. Around z ~ 200, the neutral gas temperature decoupled from the 
CMB and adiabatically cooled in an expanding Universe. During this period (lasting to z ~ 30), 
called the Dark Ages, the 21-cm line can be seen in absorption and the physics of this period is 
relatively well-understood. At the end of the Dark Ages, the spin-temperature coupling to the CMB 
started to dominate over its coupling to the cold gas, increasing the spin-temperature such that its 
differential brightness temperature approached zero once more. Around the same time, at z ~ 30, 
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Figure 1: A graphical overview of the main phases of neutral and ionized hydrogen in the Universe. From 
Robertson et al. (2010). 

haloes of ~ 10 6 solar masses started to form in sufficiently large numbers. Gas could cool suf¬ 
ficiently in to their potential wells to form the first (Pop-Ill) stars. These stars both radiated and 
heated their surrounding gas. Via the Wouthuysen-Field (W-F) effect, the spin-temperature was 
able to couple to the cold gas temperature, creating a second period (the Cosmic Dawn ) where the 
21-cm line was seen in absoiption. While star-formation continued, it is thought that the first X-ray 
emitting sources appeared, heating the neutral gas and thereby (again via the Wouthuysen-Field 
effect) raising the overall spin-temperature to a level at or above the CMB temperature. Conse¬ 
quently, the 21-cm line became visible in emission around z ~ 15. While radiating and heating, 
these first sources (possibly including mini-quasars) also ionized the gas around them and a pe¬ 
riod of Reionization started that is thought to have lasted until z ~ 5 — 6. During the latter phase 
most of the neutral hydrogen became ionized and is thought to still reside in the IGM today. Dur¬ 
ing the same period the gas was metal enriched and seed black-holes grew to super-massive black 
holes already seen at redshifts as high as z ~ 7. As discussed in the accompanying chapters in 
this volume, observing the Dark Ages, Cosmic Dawn and Epoch of Reionization is crucial for our 
understanding of the Universe which we see from the present day up to high redshifts and provide 
an enormous potential for synergy with other IR to sub-mm facilities (e.g. Euclid, ELT/TMT/GMT, 
JWST, ALMA, etc). As such, their study features high on nearly all future science and instrumental 
roadmaps, including that of the SKA, and is one of the main science drivers for a range of SKA 
precursors and pathfinders (e.g. LOFAR, MWA, PAPER and GMRT, and the planned HERA 1 ). 

2.2 Current constraints on the Cosmic Dawn and Epoch of Reionization 

Even though the history of neutral hydrogen is crucial to our understanding of the high-redshift 
Universe, very little is known about it. There is only a handful of observations that currently 
constrains physical models of the evolution of neutral hydrogen: 

• the Gunn-Peterson trough at/near the end of reionization from which (a limit on) the HI 
optical depth can be inferred (e.g. Becker et ah, 2001; Fan et ah, 2002, 2006), 

• the discovery of Gamma-Ray Bursts (GRB) during and after reionization (e.g. Ruiz-Velasco 
et ah, 2007; Sparre et ah, 2014) from which the high-mass star-formation rate (SFR) can be 
inferred, 

1 http://reionization.org/ 
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• the metal abundances and old stellar populations (e.g. Simcoe et al., 2012) possibly in Damped 
Lyman Absorbers (DLAs), 

• the evolution of the IGM temperature after reionization from which a nominal reionization 
redshift can be inferred, if the IGM cooled adiabatically and did not re-heat (e.g. Theuns 
et al., 2002; Bolton et al., 2010), 

• the polarization of the CMB due to ionized gas (basically electrons) integrated along the 
line-of-sight to the present day from which a median reionization redshift and duration can 
be inferred (see e.g. Robertson et al., 2015), 

• the integrated IR and X-ray backgrounds setting limits on the emission from the first (Pop¬ 
ulation III & II) stars, X-ray sources and quasars (e.g. Zemcov et al., 2014; Fernandez & 
Zaroubi, 2013; Helgason et al., 2014) 

• the observations of the first galaxies through drop-out techniques, currently out to z ~ 12 
already (e.g. Ellis et al., 2013), or Lyman-a emission (e.g. Matthee et al., 2014), providing 
limits on ionizing (stellar) sources, and the HI content in the IGM from Lyman-a emitters. 

Although theoretical models are highly degenerate especially at increasingly larger redshifts, it does 
appeal - from these observations that reionization is roughly half-way at z ~ 10, followed by a rapid 
increase in the SFR-density below this redshift (e.g. Lorenzoni et al., 2011; Bouwens et al., 2012; 
Ellis et ah, 2013). Around z ~ 5 — 7 potentially 10% of hydrogen could still be in neutral patches 
(e.g. Mesinger, 2010), although most of it was (re)ionized by then, making most of the the universe 
transparent to uv-radiation all the way to the present day. Beyond z ~ 10 even less is known and 
most of what we think was happening is actually based on theory (e.g. Pop-Ill stars) and simulations 
that are constrained by physical models extrapolated from lower redshifts (e.g. X-ray heating via 
XRBs; Fialkov et al., 2014b). It it therefore crucial that new observational avenues are explored to 
open up these early phases in the Universe for observational studies. The redshifted 21-cm emission 
of neutral hydrogen provides just such a possibility. Neutral hydrogen is all-pervasive until its 
(re)ionization completes and it complements (i.e. anti-correlates) with sources of ionization, with 
the CMB and with the IR/X-ray background (although see Fialkov et al., 2015). Moreover due to 
the line-nature of the 21-cm hyperfine transition, redshift (and velocity) information can be gained 
directly from 21-cm emission by fine-tuning observations to different radio frequencies [i.e. to 
1420/(1 + z) MHz], Hence tomography of (redshifted) 21-cm emission can be done by covering a 
wide range of frequencies. 

In the case of the SKA, frequencies of 50 MHz (during phase 1) and higher allow HI mapping 
at redshifts below z ~ 27 (~ 120 Myr after the Big Bang) all the way to the end of reionization 
z ~ 5 — 6 (~ 1 Gyr after the Big Bang). Especially the highest redshifts are extremely hard to reach 
with other instruments, even with the James Webb Space Telescope (JWST) in the next decade. 

2.3 Current 21-cm detection experiments 

There are two general approaches to observe HI emission via the 21-cm line: (i) its spatially- 
averaged global signal and (ii) its spatial fluctuations, both as function of frequency (i.e. redshift. 
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time, distance). In the former case the average redshifted 21-cm emission 2 at a given redshift 
is measured over a very large (possibly global) area of the sky such that spatial fluctuations arc 
averaged away, in general via single dipoles or fully-filled antennae (e.g. dish, aperture arrays with 
near unity filling factor). The global signal can vary between a few mK brightness temperature 
near the end of reionization and around those redshift where the HI spin temperature and the CMB 
were almost same, to as much as +30 mK when the spin temperature was well above the CMB 
temperature, and to as low as — lOOmK or even less during the Cosmic Dawn when HI was seen 
in absoiption against the CMB and the spin-temperature coupled to the cold neutral gas. Whereas 
in principle these signals could be detected in a few hours to a few days based purely on thermal 
noise limits, such a detection turns out to be extremely difficult to the high dynamic range that 
is needed in band-pass calibration (i.e. > 10 6 ), in the presence of RFI and the complex spatially- 
temporal-frequency varying sky that couples to a generally polarized and frequency-dependent 
receiver beam. Whereas many global experiments arc ongoing (e.g. EDGES, CORES, SARAS, 
LEDA, etc; see Subrahmanyan et al. 2015) none have reached the few to tens of mK level yet 
required for a solid detection although bandpass stability of ~ 1 K has now been reached (Patra 
et ah, 2013, 2015). 

An alternative/complementary approach to detect the redshifted 21-cm emission is via an inter¬ 
ferometer, which is insensitive to the global signal (which is detectable on the (near)zero spacings 
only), but can only detect fluctuations in the 21-cm brightness temperature. Although these fluctu¬ 
ations can be considerable, they arc far smaller in amplitude than the global signal hence requiring 
long ( 1000 hrs) of integration time to detect. At lower redshifts during the Epoch of Reion- 
ization, however, ionized bubbles could be —30mK in depth and detectable on somewhat longer 
baselines with aremin resolution. Apart from the ionized bubbles, brightness temperature fluctua¬ 
tions arc sourced by density fluctuations in the time-evolving dark-matter distribution and via the 
spin-temperature coupling to the neutral gas and/or the CMB, though the Wouthuysen-Field effect 
and through heating. In addition, fluctuations arc sourced by peculiar and bulk-flow velocities that 
Doppler-shift the 21-cm line, affecting its brightness temperature. Whereas this complex astro¬ 
physics make interpretation of the 21-cm signal harder, it also provides a wealth of information on 
sources and physical process occurring during that period. In combination with the global signal 
and other observables (e.g. galaxies, CMB polarization, etc) these brightness-temperature fluctu¬ 
ations provide a treasure trove of information. The second advantage of using an interferometer 
is to reduce the levels of foreground emission from the Milky Way (several Kelvin versus several 
hundred Kelvin) that could potentially contaminate the feeble 21-cm signal, as well as the fact that 
interferometers are often easier to calibrate than total-power experiments. At the moment each 
of these approaches arc pursued by multiple teams and valuable in their own right and targeting 
different parameter spaces, thereby complementing each other. Current HI detection experiments 
arc ongoing with the GMRT (e.g. Paciga et al., 2011), LOFAR (e.g. Yatawatta et al., 2013; van 
Haarlem et al., 2013), PAPER (e.g. Ali et al., 2015; Pober et al., 2015), MWA (e.g. Dillon et al., 
2014) and observations with the LWA (i.e. LEDA) (e.g. Greenhill et al., 2014), as well as with 
NenuFAR 3 , under construction, are planned. 


2 More precisely its brightness temperature. 

3 http://nenufar.obs-nancay.fr/ Argumentaire-scientifique.html 
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2.4 Planned high-redshift 21-cm arrays 

Despite great strides forward over the last decade, no high-redshift 21-cm detection has yet been 
successfully claimed, neither of the global HI signal nor of its fluctuations. Even if successful, all 
current observational programs (except for those targeting the global signal) aim for a statistical 
detection via power-spectra (and high-order statistics) and only on the very largest angular scales 
(degrees) could one potentially reach a S/N~few in a few MHz frequency bin, comparable to 
the first CMB maps with COBE (Smoot et ah, 1992). To overcome the current S/N~l barrier 
substantially more collecting area is required, especially on short baselines corresponding to a few 
to tens of arc minute scales, by an order of magnitude over the current LOFAR-core collecting area, 
the latter being the largest low frequency array currently aiming to detect high-z 21-cm emission. 
Such a collecting area is foreseen in the current baseline design of SKA1 and might evolve in to 
SKA2 with four times the SKA1 collecting area (see this volume), allowing one to move from 
upcoming statistical detections to tomographic (i.e. direct imaging) measurements of HI in the 
next decade. A similar US-lead effort called HERA 4 is also planned extending the current PAPER 
project in South-Africa, and NenuFAR in France that aims for a nearly fully filled low-frequency 
(10-80MHz) array of 400-m diameter. 

3. Some relevant physics during the CD/EoR eras 

Radio-telescopes measure the sky’s intensity distribution (more precisely its Fourier transform), 
which traditionally is expressed in terms of a brightness temperature T b in the Rayleigh-Jeans 
regime. It can be shown that the brightness temperature of neutral hydrogen as function of redshift 
z, seen against the CMB, can be written as (e.g. Madau et al., 1997): 


8T b 
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1 +z — 7cmb \ 

"hT V T a J 


(AUU 

V 1-0.248 J 



4 £2 b h \ 
v 0.044 0.7 J 
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H(z) dru) 



l 


mK. 


(3.1) 


This equation contains several terms set by cosmology (i.e. the global baryonic density and 
spatial density fluctuations, fU and 5, respectively; the total mass density and the Hubble con¬ 
stant £2 m and H(z) or h , respectively), by (g)astro-physics (i.e. the spin-temperature T s of HI, 
the CMB temperature 7cmb) and by Doppler effects (dvu/drn; e.g. due to peculiar motions and 
bulk-flows of the gas). Via the brightness temperature and different measures of it, e.g. variance, 
higher-order statistics, power-spectra, n-point correlations, tomographic cubes (i.e. images), HI- 
absorption spectra, cross-correlation, etc., one can address questions about the physics and sources 
responsible for processes during the CD/EoR such as heating (changing T gas ), Lyman-a emission 
and the Wouthuysen-Field effect (changing T s via r gas ; Wouthuysen (1952); Field (1958, 1959)), 
reionization (changing xhi) and the growth of density fluctuations and peculiar velocities (changing 
8 and dvu/dr^). 

4 http://reionization.org/papers/ 
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3.1 Phases of 21-cm emission and absorption 

Largely due to the expansion of the universe and within it the gravitational collapse of matter, the 
universe at some point becomes non-linear in its (evolving) density structures producing the first 
stars, stellar remnants, black holes, and galaxies. These objects in turn will impact the IGM/HI 
surrounding them, leading to a host of observational effects that can be used to assess their proper¬ 
ties. Below we shortly summarize the main phases of HI after recombination (z ~ 1100) and before 
complete reionization (z ~ 5 — 6). We again refer to Barkana & Loeb (2001); Furlanetto & Briggs 
(2004); Lewis & Challinor (2007); Pritchard & Loeb (2008); Morales & Wyithe (2010); Loeb & 
Furlanetto (2013); Fialkov et al. (2014a) for many more details. 

• Age of Ignorance (z~1100-200): Just after recombination the spin-temperature of neutral 
hydrogen tightly coupled to the CMB temperature, indirectly via the gas temperature (via 
Compton scattering), making the 21-cm line invisible against the CMB temperature. At the 
moment no known observational technique is able to probe this era, even if it were possi¬ 
ble at these high redshifts. Observing these high redshift using the 21-cm line would first 
of all have to be carried out in space to avoid the ionospheric plasma frequency, but most 
likely would be limited by self-absorption in the MW (Novaco & Brown, 1978). In case 
of dark-matter annihilation (see e.g. Valdes et al., 2013) or other heating sources (Lyman-a 
photons from recombination; (see Fialkov & Loeb, 2013)) could one conceivably have the 
spin-temperature decouple from the CMB temperature. However, it is likely that this period 
in the history of the Universe will remain out of reach till far in the future, and remain a true 
Age of Ignorance. 

• Dark Ages (z~200-30): About five million years after recombination the spin-temperature 
of hydrogen coupled - via trace electrons - to the colder-than-the-CMB gas-temperature. 
Neutral hydrogen became visible in absorption again the CMB during that era. At the end 
of the Dark Ages the density of trace electrons, however, dropped sufficiently (through ex¬ 
pansion of the Universe) to make the coupling of the spin and gas temperature inefficient. 
The spin temperature started to follow the CMB temperature again and 7j, approached zero 
from below. However, this phase just preceded the formation of the first radiating sources 
and probably only lasted briefly or might not even have been fully reached (i.e. 7j, ~ 0 mK). 
Processes that can be studied through measurements of 7], during the Dark Ages are the dark- 
matter power-spectrum evolution and its annihilation physics, baryonic bulk-flows (Tseli- 
akhovich & Hirata, 2010) and the physics of gravity and general relativity. The physics 
during this era can largely be understood through linear theory (Lewis & Challinor, 2007), or 
lineal - corrections thereof (e.g. Ali-Halmoud et al., 2014) and deviations of observation from 
CDM predictions will immediately indicate new physics. Studying HI deep in to the Dark 
Ages however will remain out of reach in the near future and requires space-based radio tele¬ 
scopes because 21-cm emission will have been redshifted to near or below the ionospheric 
plasma frequency. 

• Cosmic Dawn (z~30-15): The Cosmic Dawn is typically defined as the time when the 
first stars (or other radiating sources) were formed. Lyman-a emission from these sources 
efficiently coupled the spin temperature to that of the cold gas (via the Wouthuysen-Field 
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effect), again leading to neutral hydrogen seen in absoiption. At the same time, however, the 
gas itself was heated supposedly via X-ray heating or, possibly, via the hot ISM and Lyman-a 
photons (Pacucci et al. 2014). This gas-heating and the parallel process of coupling the spin 
and gas temperature led to a rapid rise in the brightness temperature of neutral hydrogen until 
it is finally seen in emission around z ~ 15. One should note that many of these processes arc 
still ill-understood and all these effects (heating, coupling, etc) could shift around in redshift 
substantially, especially between the Cosmic Dawn and the Epoch of Reionization. Hence 
when these processes exactly occurred is not known and redshifts indicated here arc merely 
indicative numbers currently expected from nominal models. Processes that can be studied 
during the CD are the formation of the first (pop-111) stars, the first BHs, X-ray heating 
sources, W-F coupling, bulk-flows, etc. 

• Epoch of Reionization (z~15-6): While heating together with the W-F effect change the 
spin-temperature of the neutral hydrogen, the same radiation field (i.e. uv-radiation) starts to 
ionize hydrogen leading to the percolation of bubbles around the first mini-haloes containing 
(Pop-Ill) stars and possible intermediate mass black holes, i.e. mini-quasars. As time pro¬ 
gresses and the universe becomes increasingly more non-linear (on small scales), more stars 
and quasars arc formed. Although recombination can have some impact (see e.g. Sobacchi 
& Mesinger, 2014), it can not stop or balance ionization in an expanding universe and by 

6 the entire universe, apart from pockets of neutral hydrogen (mostly in galaxies), will 
be ionized once more. Processes that can be studied during the Epoch of Reionization are 
the physics of the ionizing sources, such as pop-111 and II stars, mini-quasars, feedback to 
the IGM and the transition to the currently visible universe. 

• Post-Reionization (z<6): Whereas most neutral hydrogen will have been ionized, pockets of 
neutral hydrogen might remain even at lower redshifts and in galaxies (e.g. Mesinger, 2010), 
which could be studied though intensity mapping and via cross-correlations with other e.g. 
IR surveys (see this volume). 

3.2 Sources of heating, radiation, feedback and ionization 

Whereas the above situation delineates general eras during which certain physical processes (e.g. 
Lyman-a emission, heating, ionization, etc) might dominate, the understanding of these processes 
is far from certain and depends strongly on the types of sources responsible for uv- and X-ray 
emission. Via their imprint on the brightness temperature of neutral hydrogen one can gain insight 
in to these first radiating sources, which is one of the main reasons to study high-z 21-cm emis¬ 
sion. Below a range of sources is listed that arc currently thought to play a role. We refer to the 
accompanying chapter in this volume for more details. 

• Population-Ill and -II Stars 

The first stars presumably formed just after the Dark Ages from relatively pristine gas (pri¬ 
mordial abundances), and will soon be detectable by upcoming missions, such as JWST, 
through their supernova explosions (de Souza et al. 2013, 2014). These first stars (Pop- 
III) are thought to have been of high mass, although recent work suggests they fragmented 
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possible into stars of several tens of solar masses (e.g. Stacy & Bromm, 2013). It is the uv- 
radiation of these stars that is thought to couple the hydrogen spin temperature to that of the 
cold gas. Emission from resulting X-ray binaries (XRB) could subsequently lead to X-ray 
heating. In addition, the first and second generation (Pop-II) stars arc thought to be responsi¬ 
ble for the ultimate (re)ionization of nearly all neutral hydrogen during the Epoch of Reion¬ 
ization. What the relative roles of Pop-Ill and Pop-II stars exactly are in these processes, i.e. 
W-F coupling, heating and ionization, is currently ill-understood and the processes involved 
in the formation of Pop-Ill stars and whether they form the seeds for SMBHs is also not 
clear. At lower redshifts where the first galaxies can be observed, however, it is clear that 
stars responsible for reionization must be in mini-haloes and in galaxies ~500 times fainter 
(~7 mag.) than current observational limits (e.g. Dave et al., 2006). Hence enormous extrap¬ 
olations arc needed to make any inference from current observations, let alone to redshifts 
even higher. Redshifted 21-cm emission observations clearly arc important in setting limits 
on star formation and the types of stars (or AGN) that partake in the early phases of the 
CD/EoR eras. 

• Mini-Quasars and AGN 

A second source of ionization and heating can be intermediate-mass black holes (IMBHs) 
and resulting mini-quasars due to accretion discs or Bondi accretion (e.g. Dijkstra, 2006). 
Whereas it is thought that mini-quasars arc not the dominant source of ionizing photons, 
they can still play a secondary role and possibly be a source of harder X-ray photons that 
can more uniformly heat the IGM well above the CMB temperature. These IMBHs are 
also thought to be the seed-BHs of present-day SMBHs and AGN, although the process of 
accretion appeal's to be super-Eddington to have them grow from 100 solar masses (expected 
from massive Pop-Ill stars) to 10 9 solar masses in the AGN seen already at z ~ 7 (Mortlock 
etal.,2011). 

• X-ray Binaries 

Great uncertainty also remains whether X-ray binaries are sources of heating at high redshifts 
and whether heating takes place on global scales. If heating is inefficient some IGM patches 
could remain very cold and below the CMB temperature causing the redshifted 21-cm bright¬ 
ness temperature to remain in absorption substantially impacting the level of strength of the 
total-intensity and fluctuation signals from neutral hydrogen. This could lead to considerable 
effects even at very low redshifts during the Epoch of Reionization (e.g. Pritchard & Loeb, 
2008). 

• Dark-Matter Annihilation 

Literately a "Dark Horse" 5 is whether dark matter could be a source of ionization or heating 
at very high redshifts during the Dark Ages, already influencing the evolution of the IGM 
well before the Cosmic Dawn (e.g Taoso et al., 2008; Ripamonti et al., 2010). 

Having outlined a number of important eras as well as the physics and source responsible for 
the brightness temperature evolution of neutral hydrogen (both spatially and in time), we now 

■’Benjamin Disraeli, in "The Young Duke" (1831) 
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continue to discuss which observables the SKA can obtain and use to constrain these processes and 
from it learn about the sources responsible. We also shortly mention possible synergies with other 
(planned) facilities. 

4. Observables via the redshifted 21-cm brightness temperature 

Redshifted 21-cm emission from neutral hydrogen manifests itself in multiple ways and various 
methods and techniques can be used to extract information from it. In the following section we 
discuss a number of such approaches, whereas in the subsequent section a short summary is given 
of other observables (and facilities) that the 21-cm brightness temperature could be cross-correlated 
with. We refer to Mellema et al. (2013), as well as to Chang et al. (2015) and Jelic et al. (2015) in 
this volume, for more details on the synergy between the SKA and other observatories. 

4.1 Measures of the brightness temperature field 

All of the information one hopes to gain from the physics occurring during the Dark Ages, Cosmic 
Dawn and Epoch of Reionization from neutral hydrogen is obtained via its redshifted 21-cm bright¬ 
ness temperature. As shown in Section |||, this temperature is set by a combination of many different 
physical processes which need to be disentangled. As of yet, it is not fully clear whether this can 
be done completely or not. In addition, it is likely that in the first experiments, also with the SKA1, 
only limited measurements of the brightness temperature will be made. Below we summarize some 
of the standard measures of 7],: 

• Total Intensity Redshifted 21-cm emission of neutral hydrogen averaged inside "shells" of 
common redshift (or observed frequency) are close to impossible to measure with interfer¬ 
ometers since they measure intensity differences (Subrahmanyan et al. 2015). However, as 
seen in Section |3], the total expected brightness temperature is of order ~ 27xhi(1 + 8) mK at 

lOif T s > 7cmb- When averaging over large areas of the sky (2> 1 degree) (8) « 0 and 
assuming no reionization (i.e. xhi = 1), one can determine the total intensity of the signal. 
This is "most easily" achieved with single receiver/dipole measurements and is discussed 
in much greater detail in Subrahmanyan et al. (2015). Measuring this signal as function of 
redshift allows one to learn about different physical processes such as Lyman-a coupling 
by the first stars, heating by (X-ray) sources, reionization, etc. At this point it is not clear 
whether the SKA will be capable of measuring this signal since it requires a very accurately 
calibrated instrument. 

• Variance or second moment The second-order effect comes from variations in the bright¬ 
ness temperature due to variations of 8, xm, E and v||. These variations lead to a spatial 
variance in E as function of observed frequency or redshift (see Mellema et al. 2015 and 
Wyithe et al. 2015) which can be measured in interferometric images, assuming that sys- 
tematics and thermal noise are under control and/or known (although thermal noise can be 
removed via cross-variance techniques over small time-frequency steps). By measuring the 
variance one can again learn how HI evolved over cosmic time, in particular when reioniza¬ 
tion peaked (the contrast between HI and bubbles is very large and leads to a large variance). 


10 



CD/EoR with SKA 


L.V.E. Koopmans 


However, the shape of PDF(7b) also teaches us about its moment (skewness and/or kurto- 
sis) again revealing information about reionization (e.g. Harker et al., 2009). In particular 
the skewness of the PDF might be a tell-tale signature of reionization that would be hard 
to mimic by systematic effects (which often are strongly correlated between frequency and 
scale smoothly with frequency as well). 


(Anisotropic-)Power-spectra One step further would be the analysis of 7j, in the Fourier 
domain in terms of the variance as function of spatial (or angular) scales, i.e. the power- 


spectrum. The power-spectrum depends on the spatially varying parameters in Equation 3.1 


being xhi, 5, T s and vy. All other parameters can be assumed constant and only enter via the 
general FRW metric. The power-spectrum is directly related to the underlying local phys¬ 
ical processes and sources of reionization (governing xhi), heating (governing T s ), as well 
as to the evolution of the density field (governing 8 and vn). To first order one can Taylor 
expand equation 3H in these parameters, Fourier transform the result, from which the power- 
spectrum can be formed. Since the parameters are not independent (e.g. an over-density can 
also lead to extra heating and ionization), the resulting 21-cm power-spectrum becomes a 
rather complex function of these parameters and their cross-correlations, and can become 
anisotropic. These can either be determined via simplified schemes or via direct numerical 
simulations (which also contain simplified sub-grid physics). Disentangling all effects from 
the power-spectrum as function of redshift is nontrivial and remains an ongoing and active 
research field. To first order, the power-spectrum is isotropic and can be spherically aver¬ 
aged (where the frequency direction over some limited bandwidth is a proxy for comoving 
distance, hence inverse wavenumber), although in the case of SKA (possibly even before) 
the sensitivity if sufficient to measure anisotropy in the power-spectrum which can be an 
exciting measure of underlying physics (e.g. Fialkov et al., 2015). Hence it is often better to 
express the power-spectrum in two dimensions, separating the wave vector in a component 
perpendicular to the line of sight and a component parallel to it, in particular to study pecu¬ 
liar velocities, the light-cone effect, as well as potentially identify still-remaining systematic 
errors in the data-processing. 


• Higher-order statistics/n-point correlations: Already touched up on above, the brightness 
temperature field after the Dark Ages will start to deviate from a Gaussian random field. This 
means that the power-spectrum or two-point correlation function will become an increasingly 
more incomplete description. The brightness temperature difference between two points will 
deviate from a Gaussian distribution, and the three and high-point statistics will become 
non-zero (Pillepich et al., 2007; Cooray et al., 2008). Studying these statistics can therefore 
provide additional information about physical processes, even if direct imaging might still 
be infeasible. 


• Tomography/Images: The ultimate goal of any HI observation is make direct images and 
determine the brightness temperature with high significance for each spatial-frequency cell 
(see Mellema et al. 2015 and Wyithe et al. 2015 this volume). This clearly contains the max¬ 
imum amount of information obtainable from the data-set, but interpretation is not always 
straightforward and the analysis of images has been a sorely neglected topic in the litera- 
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ture. In fact, even from direct numerical simulations, results are often presented in limited 
statistical measures such as the power-spectrum. Developing new ideas on how to analyze 
brightness temperature images is therefore crucial before the SKA comes online since it will 
be the first instrument capable of imaging/tomography on scales of several arc minutes to 
degrees. Only in exceptional cases might one expect to do this pre-SKA. 

• Bubbles and HI topology: Whereas initially imaging of HI brightness temperature fluctua¬ 
tions at the level of several mK might not be feasible (e.g. in roll-out, etc), the stark contrast 
of 20-30mK and size (ten(s) of Mpc) of ionized bubbles will probably allow them to be seen 
relatively easy on during the roll-out and later phases of SKA1 imaging. Their distribution 
as function of size, shape, etc. will likely reveal swiftly how reionization proceeded (i.e. 
inside-out, or otherwise). It might also reveal the nature of the sources causing reionization, 
in particular - if (mini-)quasar play a role. Although some of this information is encoded in 
the power-spectra (which forms from a combination of the non-ionized HI field times a mask 
where HI has disappeared) their precise sizes, shapes, etc. will provide a far richer source of 
information. In addition it might be possible to study the local velocity field structure around 
bubbles which might reveal additional information that the power-spectra do not. 

• 21-cm absorption: Whereas all the above measures of neutral hydrogen are based on its 
brightness temperature contrast with the CMB, there could be radio sources (mJy or brighter) 
present at sufficiently high redshift that the 21-cm line could be see in absorption against 
such source (e.g. Carilli et al. (2002) and Ciardi et al. 2015). A mJy source (e.g. AGN or 
GRB) of a few arcmin resolution at say 150MHz would already be sufficient to measure the 
HI forest and from it (equivalent to the Lyman-a forest) measure the line-of-sight power- 
spectrum of neutral hydrogen. Although success depends on the (still unknown) presence 
of radio sources at redshifts well beyond the end of reionization, it would provide a unique 
opportunity to measure /F| (k) on k-modes far larger (i.e. scales far smaller) than accessible 
via imaging/tomography or power-spectra as discussed above. As such SKA1 and 2 could 
provide the first measures of mini-haloes during the Epoch of Reionization something that 
even JWST will find difficult to do. 

4.2 Synergy with other surveys/Cross-correlations 

Whereas SKA will characterize the brightness temperature in detail via either power-spectra, direct 
imaging or 21-cm absorption possible to redshifts of z ~ 27 during the Cosmic Dawn, and will 
be unique at this, complementary observations can further enhance the capability to extract infor¬ 
mation from this rich data-set. Among these are: (1) the Planck (or future) CMB maps, (2) the 
Far-Infrared background or emitters, (3) surveys of high-z Lyman-a emitters, (4) high-z Lyman- 
break Galaxies/Dropouts, (5) high-z SNae/GRB transients, (6) CO/CII emitters and (7) DLAs. We 
refer to Chang et al. (2015) and Vibor et al. (2015) in this volume for details. 

5. Survey Design 

In this section we outline a potential observational survey with SKA1 (and 2) to characterize the 
redshifted 21-cm brightness temperature of HI between z = 27 (v=50MHz) and z~6 (end the EoR), 
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via its two main observable: power-spectra and tomography. We assume throughout that calibration 
and foreground subtraction is done to a level accurate enough that it does not affect the resulting 
residual data-cubes (containing both the 21-cm signal as well as noise). For more details on the 
latter we refer to Mellema et al. (2013). Hence thermal noise and sample variance (in the case of 
power-spectra) arc assumed to dominate the error budget in the power-spectra and thermal noise 
dominates the error-budget in the images. 


5.1 Scaling Relations for power-spectrum sensitivity 


To understand the choices of survey specifications we need some guidance. This is most easily 
understood by assuming a uniform nv-coverage in an interferometric experiment. As shown in 
Mellema et al. (2013) and based on McQuinn et al. (2006), the power-spectrum error due to thermal 
noise is then given by: 


A 


2 

Noise 



eA e ff 


A 2 

^coll 


(5.1) 


The total error is obtained by adding the sample variance to this, which depends on the power- 
spectrum of the signal itself and the number of independent samples per volume. Whereas this 
equation gives a reasonable idea of the level of thermal noise variance on 21-cm power-spectra, it 
can be tilted as function of k = Iti/L, depending on the nv-density as function of baseline length. 
In this equation, the co-moving distance to the survey redshift is D c , the depth of the survey for 
a bandwidth B is A D c and the field-of-view of the survey is Gp 0 y. In case of multi-beaming, the 
number of independent beams is A),. A system temperature of T sys is assumed and an integration 
time q nt . Each station/antenna is assumed to have an effective collecting area A e rr, that corresponds 
to the field-of-view as QpoV = A 2 /A c ff, where A is the wavelength at the center of the frequency 
band. The antennae are distributed over a core area A core in such a way that the nv-density is 
uniform and the total collecting area (the sum of all effective collecting areas of each receiver) is 
assumed to be A co u. The equation has been shown to be correct via numerical integration, using 
the method outlined in McQuinn et al. (2006). Regardless the assumption of a uniform nv-density 
(which affect the exponent in fc 3 / 2 and the normalization), the scaling with the layout of the array 
follows 


A 2 - 

^Noise 


sV^eff 


A 2 

A coll 


o^stat^eff 


3/2 


, V^stat^coll, 


(5.2) 


where we assume A eo u = /V staL Aeir with AQat stations/receivers. Given the above equation and scaling 
relations, a number of conclusions can be drawn for the power-spectrum sensitivity at a given k- 
mode scale: (i) sensitivity scales most rapidly with collecting area, (ii) increasing field-of-view 
helps, but only with the square root of the field of view, and (iii) more compact arrays (still covering 
the /c-modes of interest) arc more sensitive. All current arrays (i.e. MWA, LOFAR and PAPER) 
follow strategies that optimize (or maximize within budgetary limits) these three parameters. SKA1 
and 2 will both increase A co n by 1-2 order of magnitude over all current arrays and minimize A core 
by maximizing the filling factor and placing much of the collecting area in a rather small core area. 
Finally, based largely on computational limitations, the field of view of SKA is at least of order 
five degrees at 100 MHz, but during later phases (e.g. SKA2) could be increased through either 
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multi-beaming or through reducing the beam-formed number of receivers. Although in considering 
optimizing an array for CD/EoR science, one should not forget that the system also needs to be 
calibrated for instrumental and ionospheric errors, and foregrounds (compact and extended, in all 
Stokes parameters) need to be removed, which might require the use of long baselines. Any array, 
also SKA1 and later 2, therefore preferably would consist of a compact core with "arms" extending 
to many tens of kilometer. 

5.2 Scaling relations for tomography/imaging 

Si mi lar to the power-spectrum sensitivity, one can also obtain a scaling relation for imaging or 
tomography. One readily finds that the thermal noise inside a resolution element corresponding to 
a scale k scales as: 



assuming a sinc-function as model for the sky-components (i.e. a top-hat in uv-space). In this case 
A e ff drops out of the equation, although we leave it in to retain a si mi lar form to the equation for the 
power-spectrum. We also see that on larger angular scales (i.e. smaller k±), the noise decreases and 
that a more compact array helps to reduce the thermal noise per resolution element, by increasing 
the number density of visibilities per nv-cell. The latter form of the equation expresses this through 
the tilling factor (j],\\) of the core area in terms of collecting area and the field-of-view in a reso¬ 
lution element of the core. Filtering longer baselines (or smoothing of the image) to a resolution 
corresponding to k averages down the noise and lowers the brightness temperature error, hence k 
does not necessarily correspond to the largest baseline afforded by the core area. The latter simply 
sets the nv-density. For example for the inner paid of SKA1 the core of 400m diameter has /mi ~ 1 
and on scales k± = 0.1 cMpc 1 (~ 20') one expect an error of ~ 0.2 mK per MHz in 1000 hrs 
of integration time at 150MHz, assuming r sys = 400 K, which is below the expected brightness 
temperature fluctuations on that scale. However, with increasing resolution the filling factor of the 
core decreases and k± increases, both leading to a rapidly increasing brightness temperature error 
ultimately limiting the angular scale where imaging can still be done. For SKA1 and 2 this limits 
occurs on an angular scale of a few to about ten arc minutes in a reasonable (~1000 hi - ) integration 
time. 

5.3 Current SKA1-LOW Baseline Design and its Relevant Design Specifications 

In the following we assume the current SKA1 baselines design (BED), based on the station layout 
provided in Braun (2014). Although the latter is not necessarily the ultimate layout, it represents 
a "close-enough" description of the array that no major differences are expected in the resulting 
sensitivity calculations, except if the array is (i) substantially thinned out (i.e. reducing the total 
collecting area), (ii) dramatic changes in antenna layout are introduced (i.e. changing the filling 
factor) or (iii) substantially reducing the instantaneous field of view (i.e. through beam forming 
larger numbers of receiver elements). We refer to other chapters in this volume for a full description 
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of the current BLD. Using these antenna positions (or radial distribution) we calculate full 12- 
hr wv-tracks in the direction of the South Galactic pole and average the nv-density within radial 
bins. Similar densities are obtained for tracks at lower elevations. We assume a nominal 35m 
beam-formed "station" as well, which leads to a frequency-dependent field-of-view, being around 
5 degrees FWHM at 100 MHz. 

The fiducial layout that we assume in this review is listed below. We emphasize that this is not 
necessarily the most optimal observing strategy, although in some cases (e.g. freq. coverage) 
full use of the capabilities of SKA1 and 2 are made. As we will discuss below, multi-beaming 
can dramatically improve the return of a single CD/EoR sutyey. Also integration times can be 
increased 6 . Finally we assume for the nominal survey that phase-tracking of a single field is done, 


rather than drift-scans although we come back to that in this Section and Section 6.3 


1. Antenna layout: We assume for SKA1 the layout as (generically) specified the BLD, as¬ 
suming that for SKA2 the collecting area is increased by a factor of four. The latter is 
implemented here via a scaling of all baselines by a factor of two and quadrupling the num¬ 
ber of stations at each SKA1 station position. This ensures that the filling factor in the core 
remains of order unity and decreases in a scale-invariant manner to SKA1. Other choices are 
possible and will be investigated in the future. 

2. Frequency range: Based on the current baseline design of SKA1, we assume a full instanta¬ 
neous frequency coverage of 50-350 MHz (i.e. z= 27-3). This range is thought to encompass 
the full Cosmic Dawn and Epoch of Reionizadon (see Mellema et ah, 2013). We assume 
that during SKA1 the bandwidth can at least be split in two beams of 150-MHz, increasing 
the survey speed by a factor two. We note again that more beams can reduce on-sky time 
significantly. In the future one might consider extending SKA1 to lower frequencies (below 
50 MHz) during the upgrade/extension to SKA2, although this will not be considered in this 
chapter. 

3. Survey FoV: The field of view of a single primary beam is assumed to be FWHMw 5° at 
100 MHz 7 . At a redshift of 

Z 10, the comoving distance is ~ 9700 cMpc, hence the co¬ 
moving scales covered by the primary beam extend up to ~1 cGpc, about ten times larger 
than the scale where the dark-matter power-spectrum peaks (~ lOOcMpc). For a bandwidth 
of 10 MHz (see below), the volume depth is ~ 170cMpc. This means that sample variance 
on scales of a few hundred cMpc is ~ 10% and less on smaller scales. On larger scales, sam¬ 
ple variance rapidly increases and only a larger survey volume (i.e. increased field-of-view 
via multi-beaming or drift-scanning) can further reduce this source of variance. One might 
consider increasing the FoV by beam forming a smaller number of log-periodic elements. 
We discuss the effect of multi-beaming/increased FoV below. 

6 We note that largely only winter-nights can be used during which the Sun and Milky Way center are below the 
horizon and the ionospheric effect are milder. Hence for deep pointed integrations, accumulating a thousand hours of on- 
sky time might take multiple years as experience with current pathfinder/precursors of the SKA already shows. Leeway 
to offset collecting area against integration time might therefore not save costs since it can lead to extremely long project 
lifetimes. 

7 We assume a D=35m beam-formed primary beam at 100 MHz, although not necessarily physical stations of 35m. 
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Figure 2: Deep Survey: Shown are the expected brightness temperature power-spectra from Mesinger et al. 
(2014) and the expected thermal-noise and cosmic-variance errors (following McQuinn et al., 2006) for the 
fiducial deep-survey parameters outlined in the text for SKA 1-LOW. The thick/thin curves are for 5 versus 
1 beam. The middle-column panels (red curves) indicate the expected total S/N ratio in such a deep survey, 
reaching a peak S/N~ 70, in good agreement with the results in Mesinger et al. (2014). The right column 
panels show the brightness temperature sensitivity as function of resolution, where the A-scale is transformed 
directly to angular scale at the corresponding redshift. The latter does not depend on the number of beams. 
The dashed line is for a BW that matches the angular scale rather than being fixed to 10MHz. The blue 
dot-dashed line is the expected variance on that scale by integrating over the power-spectrum. 
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4. Integration times: In this review, we assume an integration time of 10, 100 and 1000 hrs, 
for the different levels in a three-tiered survey with decreasing survey area, respectively. The 
reason for this choice is that SKA1 can reach levels of ~1 mK in 1000-hrs integration time 
(per 1-MHz bandwidth) over a wide range of angular scales. However, this integration should 
only be thought of as a reference, since the level of brightness temperature can vary between 
models. This integration time is also often used in the literature as a reference, making 
comparison between models easier. The overall on-sky time of the assumed (three-tiered) 
survey, however, will be larger (see Section|6j). 

5. Frequency Bandwidth/Redshift Binning: We assume a bandwidth of 10 MHz, per red- 
shift bin (not in total!), to maximize sensitivity but minimize light-cone effects, i.e. the im¬ 
pact of the evolution of the hydrogen brightness temperature with redshift mixed with the 
assumption that redshift can be regarded as a third spatial dimension of the data volume. The 
assumed bandwidth is a reasonable compromise. We note that bandwidth enters both in the 
survey depth (via AD c ) and in the sensitivity per visibility. 

6. Multi-Beaming: We assume a single primary beam (N\, = 1) for a 35m station in SKA1 for 
the full bandwidth (BW), possible two when split in to two 150-MHz bands, increasing for 
SKA2 potentially to covering the full log-period dipole beam (Ab +> 1). The latter would 
require an increase in correlator capacity beyond the currently envisioned SKA1 correlator. 
Below we will show the impact of multi-beaming, in particular in dramatically reducing 
sample variance on large angular scales. 

7. System Temperature: We assume T sys = 100 + 400 x (v/150MHz) -2 - 55 K. Although this 
can change from field to field and can increase/decrease towards/away from the Galactic 
plane, we assume it to hold to reasonable order and is typical for assumptions in the literature. 

8. Inner-Cut uv-plane: Finally we assume a 30-A inner cut to the uv -plane (i.e. the sky- 
projected baselines in units of wavelength) to reduce the impact for foreground leakage in to 
the power-spectrum. This starts limiting scales k < 0.02 cMpc 1 in the power-spectrum (See 
Fig.§) corresponding to scales of about 2 degrees or ~ 300cMpc at 150MHz. Since these 
scales already well exceed the peak of the DM and brightness temperature power spectra, 
they provide some information about cosmology, but far less information of the physics of 
CD/EoR. This cut is not needed if bright and large scale foregrounds can be controlled. 

Whereas these assumptions, based on the BLD of SKA1 and 2 should not be regarded as definitive, 
they are a guide to typical observations specification, also currently used e.g. by MWA, LOFAR 
and PAPER, and thus serve as a reference. 

6. A Three-tiered Cosmic Dawn and Epoch of Reionization Survey 

Having specified our assumed nominal observational settings, in this section, we introduce an out¬ 
line for a three-tiered survey of redshifted 21-cm emission during the CD/EoR that can address the 
wide range of physics questions posed in Section |3| 
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6.1 Main Observation Targets: 


It should be stressed that observations of HI emission from the CD/EoR have thus far not been 
made (although upper limits are currently being set) and any present-day expectations are based 
on relatively sparse observations in different wavelength regimes or through very different observ- 
ables(see Section 2.2 ). Any CD/EoR-Hl observational strategy should therefore cover as much of 
parameter space as possible (e.g. redshift, angular and frequency scales, field of view, etc), pro¬ 
viding for surprises that the Universe might put up on us, obviously bounded by any reasonable 
estimate of what is possible. Overall there are three general approaches to detecting and character¬ 
izing neutral hydrogen: 


• via direct imaging (i.e. tomography) of neutral hydrogen down to the 1 a « 1 mK brightness 
temperature level on >5-arcmin scales at z ~ 6, rising to degree scales at z=28. 

• via statistical (including power-spectrum) methods to variance levels la ss 0.01 — 1 mK 2 
over the redshift regime z=6-28, respectively, at k<0.1 Mpc 1 and to k<l Mpc 1 at z < 15. 

• via 21-cm absorption line observations against high-z radio sources, if present, with 1- 
5 kHz spectral resolution (2-10 km/s at z=10) with S/N>5 on 1 mJy. 

• Although SKA2 can attain SKAl’s limits 4x faster or, in the same integration time go 16x 
deeper (for power-spectra), new discoveries with SKA1 are expected to guide SKA2 obser¬ 
vations, leading to new observational targets. SKA2 can also image smaller ionized bubbles 
and sub-mK brightness-temperature fluctuations during the Cosmic Dawn, unfeasible with 
SKA1. An extension of the frequency range from SKA1 to SKA2 could also enable one to 
probe the late Dark Ages at z ~30-40. 


The general target can be reached in a set of surveys and target fields in three tiers (depth versus 
area), which we will discuss in the following three sub-sections. 

6.2 Deep Survey 

The direct science goals of a Deep CD/EoR survey are to detect and characterize ionized structures 
and HI brightness temperature fluctuations on 5-300 arcmin scales (varying) over the EoR/CD red¬ 
shift range z = 6 — 28 to 1-mK brightness temperature level and from it derive the state, thermal 
history and chemistry of the IGM, study the first stars, black holes and galaxies and constrain cos¬ 
mology, the physics of dark matter and gravity. To reach this goal, deep lOOOhr integrations are 
considered on 5 separate 20-square degree windows covering a total of 100 square-degree on-sky 
area, using the 50-200MHz frequency range with 0.1 MHz spectral resolution (for science), using 
the full SKA1-LOW array, using multi-beaming with /V h = 2 lowering the on-sky time from 5000 to 
2500 hrs. This will not be feasible with any current (or funded) array with S/N>1 and is unique to 
SKA1-LOW. This survey lowers the (expected) thermal noise ~ 10 x over LOFAR and allows di¬ 
rect imaging with S/N3> I. The Cosmic Dawn remains inaccessible to current/funded instruments, 
until this deep survey with SKA 1-LOW. The sensitivities for power-spectra and imaging are shown 
in Fig. reffig:powerspectral. 
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6.3 Medium+Shallow Survey 

The direct science goals of a Medium+Shallow CD/EoR survey are to detect and characterize 
the 21-cm power-spectrum with a peak S/N~70 measured over k = 0.02 — LOOM pc 1 over the 
CD/EoR redshift range z ~ 6 — 28 and from it derive the state, thermal history and chemistry of the 
IGM, study the first stars, black holes and galaxies and constrain cosmology, the physics of dark 
matter and gravity. A medium-deep pointed survey of 50 times 100-hr integrations would cover 
1,000 square degrees, and a preceding shallow pointed/drift-scan survey of 500 times 10-hr inte¬ 
grations would cover 10,000 square degrees, using the 50-200 MHz frequency range with 0.1 MHz 
spectral resolution (for science) and using the full SKA1-LOW array. Both survey can be carried 
out in 2500 hours each, assuming two beams of 150 MHz. Again multi-beaming could reduce 
the on-sky time substantially. A gain in power-spectrum sensitivity would be 1-2 orders of magni¬ 
tude over current arrays (i.e. MWA, PAPER, LOFAR) in the sample-variance limited S/N regime 
(, k < 0.2 Mpc -1 ), but scales k = 0.2 — 1 Mpc~* are likely to remain inaccessible until SKA1-LOW 
is build, because of severe thermal-noise limitations. Again the Cosmic Dawn remains inaccessi¬ 
ble to current/funded instruments, until SKA1-LOW is realized. The sensitivities for power-spectra 
and imaging are shown in Fig.||. 

6.4 Absorption at 21-cm 

The science goal is to obtain 21-cm absoiption line spectra from (rare) radio sources at z > 6 
to probe very small scales (i.e. k ~ 1000 Mpc 1 ) or virialized structures (mini-haloes) and from 
it derive the state, thermal history and chemistry of the IGM, study the first stars, black holes and 
galaxies and constrain cosmology, the physics of dark matter and gravity. This can be done through 
deep 1000-hrs integrations on selected sources with a flux of > 1 mJy (possibly in the Deep Survey 
fields; see Ciardi et al. 2015) with a spectral resolution of a few km/s (1-5 khz at 150 MHz) over 
the 50-200 MHz frequency range (z ~ 6 — 28). We note also here, that no observations of this kind 
have ever been done and are not expected to be feasible before SKA1-LOW is build, because of 
the same severe thermal noise limitations that currently limit direct imaging. 

6.5 Deep to Shallow Surveys: Multi-beaming and/or Drift-Scans 

Optimal observing strategies will depend on many variables: ionospheric conditions, minimization 
of side-lobe leakage of strong sources, field in the zenith, low galactic brightness temperatures and 
polarization, optimal nv-coverage, multi-beaming versus frequency coverage, inclusion of long 
baseline information, monitoring of the station beams, RFI excision, drift-scans versus tracking, 
etc. For SKA1-LOW the field of view is typically 5-10 degrees and hence the sky drifts (see Trott, 
2014, for details) though a stationary beam (i.e. non-tracking) in about 0.5-1 hr. Deep integra¬ 
tions of 1000+ hrs would therefore require 1000-2000 nights for every single deep field, which is 
prohibitively long. Tracking the field over A hrs, would allow the same integration time to be accu¬ 
mulated in N times fewer nights. Hence it appears that drift-scans would only be useful for larger 
scale shallow surveys, where the total integration per field can be small. On the other hand, because 
SKA1-LOW allows for two beams in the baseline design, one could conceive of a combination of 
tracking-drift-scan, where the field drifts from one stationary beam to a second stationary beam and 
only every 0.5-1 hr would one switch the phase center of the trailing beam to a new phase-center. 
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Figure 3: Deep-Medium-Shallow Surveys: Shown are the expected brightness temperature power-spectra 
from Mesinger et al. (2014) and the expected thermal-noise and cosmic-variance errors (following McQuinn 
et al. 2006) for the fiducial survey parameters outlined in Section 6 for SKA1. The middle panels (thick red 
curve) indicate the expected total S/N ratio in such a survey, reaching a peak S/N~ 70, in good agreement 
with the results in Mesinger et al. (2014). From top to bottom, the red curves for the deep, medium and 
shallow surveys are shown, clearly showing an increased gain in sensitivity at the larger k modes in the 
deep survey due to low thermal noise, and an increased gain in sensitivity at the smaller k modes in the 
medium+shallows surveys due to lower sample variance, leading to a relatively flat S/N curve over one dex 
in k space. For completeness, the right column show the brightness temperature sensitivity as function of 
resolution, where /.'-scale is transformed directly to angular scale at the corresponding redshift. From bottom 
to top are the deep, medium and shallow surveys. The dashed line is for a BW that matches the angular 
scale rather than being fixed to 10MHz. The blue dot-dashed line is the expected variance on that scale by 
integrating over the power-spectrum. 
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Since these ideas have not been worked out in detail (but see Trott, 2014), the precise observing 
strategy remains an open question and will depend on a combination of nights available, sensitivity, 
resolution in the image, calibration errors, etc. 

7. Role-out of SKA1 and expansion to SKA2. 

A great advantage of radio interferometers is that it can observe even if incomplete in its number of 
receiver elements and/or stations/dishes. To enable high-impact science even with an incomplete 
array it is important to think about rollout from an incomplete SKA1 to a fully complete SKA2. 

7.1 Role-out from SKA1 to SKA2: 0.5,1, 2 and 4 times SKA1 

In this section we shortly discuss three stages: (i) roll-out of half of SKA1, (ii) SKA1 and (iii) 
rollout to a full SKA2, being four times the collecting area of SKAl. We only compare these 
stagesca assuming a scaling in the number of stations (i.e. collecting area) and no other specification 
(e.g. frequency coverage, multi-beaming, etc.). The ratio of the number of stations is assumed to 
be f s where / s = 1 for SKAl. Since neither half of SKAl nor SKA2 have any specified baseline 
designs, we have to make an assumption. Based on the scaling relations in Sections |5| a high filling 
factor is critical to reach surface brightness levels of ~ 1 mK. The simplest assumption is therefore 
to take the current SKAl BLD and scale the number of visibilities (i.e. the voltage/electric-field 
coherence between two receivers) by / s 2 , but keeping the radial uv-density profile scale-invariant. 
The above scalings ensure that the filling factor at the center of the array remains the same during 
rollout which ensure maximum sensitivity on large angular scales. In Fig.|] we show the impact of 
/ s = 0.5,1.0,4.0 on the power-spectra and brightness temperature sensitivity for tomography. We 
conclude a few things in comparison to the SKAl baseline design: 

• For k < 0.1 Mpc -1 there is no power-spectrum sensitivity loss in a single beam (~ 20 square 
degrees) because one is sample-variance limited. Above k > 0.1 Mpc -1 , however thermal 
noise dominates the error budget and one quickly gains S/N. Hence SKAl in roll-out could 
already carry out a shallow wide-held survey focussed on k < 0.1 Mpc -1 . 

• During the EoR and especially during the CD, SKAl to SKA2 start to reach k > 1 Mpc -1 
with good S/N, which is the region where the power-spectrum can rise quickly due to physics 
on small angular scale (mini-haloes) and where small ionized bubble will exhibit features in 
the power-spectrum. We note an almost order of magnitude increase in sensitivity from 
SKAl to SKA2 at /-modes beyond Z ~ 0.2 Mpc -1 

• Whereas half SKAl will enable tomography on ~ 5 arcmin scale during the EoR, during 
the CD this ability is lost and only with some luck scales larger than half a degree might be 
imaged. SKAl and in particular SKA2, however, will enable tomography during the EoR, 
CD and late Dark Ages (the latter only for SKA2). As shown in Sections ^ & |6[ sheer 
collecting area and its resulting instantaneous sensitivity is needed. This can only partly 
be compensated by integrating longer, but at the expensive of substantially increasing the 
project duration (possibly to a decade or more). 
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Figure 4: Rollout from half to four times SKA1: Idem to Fig.2 & 3. The solid red lines from bottom to 
top (middle column) show the S/N for a deep survey for / s = 0.5,1.0,4.0. In the right column the brightness 
temperature (from bottom to top) is given for SKA2, SKA1 and half SKA1. 


A remarkable result is that for most of the redshift regime (z ~ 9 — 15) cosmic variance limits the 
survey below k ~ 0.2Mpc '. We note that because both the thermal and cosmic-variance errors 
scale identically with number of beams (see text), this transition scale is invariant under increasing 
the number of beams (i.e. multi-beaming), despite that the S/N increases as \J /Vbeam- 
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8. Summary, Conclusions and Highlights 

SKAl-LOW and finally SKA2-LOW will be a transformational instrument in the study of the 
physics of the first stars, galaxies and black holes, dark-matter, the IGM, cosmology and possibly 
gravity during the first one billion years of the Universe, via statistical and direct measurements 
(e.g. power-spectra and tomography) of the redshifted 21-cm line from neutral hydrogen. 

SKA1 will enable brightness temperature levels of ~1 mK to be reached in 1000 hrs of integration 
(BW=1 MHz) from z ~ 6 to z ~ 28 (>50MHz) over increasing angular scales ranging ~5-300 ar- 
cminutes. SKA2 could potentially enable even higher redshift observations, stalling to probe the 
late Dark Ages, as well as probe a much larger area of the sky (via multi-beaming) to a much 
greater depth in less time, reducing the thermal noise variance that will still limit SKA1 on the 
small scales where the dominant source population (e.g. mini-haloes) have their largest impact 
on neutral hydrogen (both via ionization and spin-temperature effects). The impact of SKA2 is 
expected in pushing tomography (i.e. direct imaging) to higher redshifts and to greater S/N (e.g. to 
probe redshift space distortions, higher order statistics, specific regions in the cubes, etc) where the 
sensitivity of SKA1 will still be insufficient. Although the Cosmic Dawn will be opened by SKA1, 
SKA2 will allow it to be studied in much great detail similar to the impact of SKA1 on the study 
of neutral hydrogen during the Epoch of Reionization. 

We have outlined a three-tiered survey (deep and shallow/medium-deep and 21-cm absorption line 
measurements) with SKA1 that allow one to reach the goals of measuring the HI brightness tem¬ 
perature to a level of around 1 mK for tomography (a genuinely unique capability of SKA1 and 
SKA2) over a range of a few arc minutes up to a degree (depending on redshift), reduce sample 
variance and reach power-spectrum measurements up to one arc min scales. Such a three tiered sur¬ 
vey (assuming 21-cm absorption measurements can piggy-back on the deep survey) can be carried 
out in 3 x 2500 hrs 8 assuming that the full 300 MHz BW of SKA1 can be split in to two beams 
of 150 MHz. Such a survey is expected to take of order five years, assuming an efficiency of 50% 
of all night-time data (assuming 8-hr tracks) over such a period by selecting the best observing 
times and minimizing the impact of the ionosphere and emission from the Milky Way. This is 
not a project that can be carried out in a few months of full-time observations. Experience with 
the low-frequency SKA precursors and pathfinders confirms that good "observing conditions" arc 
necessary for a successful CD/EoR projects. 

Finally, a three-tiered survey with SKAl-LOW is expected to generate strong synergy with other 
major (future) facilities such as JWST, ALMA, Euclid, as well as with (CO/CII) intensity-mapping 
arrays and large ground-based IR facilities such as the ELT. Many of the details can be found in the 
accompanying CD/EoR chapters in this volume. 

8.1 SKA1 and 2 Science High-lights 

Some high-light science results, in random order, expected to be largely unique for SKA1 (and 2) 
arc given below. For each bullet point we refer to the first author of the related chapters in this 
volume: 

8 We note that collecting 2500 hrs of "good" deep pointed data might take up to ~5 years since winter-nights where 
the MW and Sun are below the horizon and the ionosphere is mild might only allow up to 500 hrs on-sky time to be 
accumulated per calendar year. Multi-beaming can significantly speed this up. 
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• Direct imaging of ionized regions and HI fluctuations on scales of arcminutes and larger 
during the Epoch of Reionization and Cosmic Dawn — Mellema et al. 2015; Wyithe et al. 
2015 

• Probing many aspects of cosmology/cosmography out to the highest possible redshifts (i.e. 
z = 27 in case of SKA1) — Pritchard et al. 2015 

• Enabling one to probe k ~ 1000 M pc 1 scales via 21-cm absorption, probing mini-halos even 
out of reach of JWST — Ciardi et al. 2015 

• Direct study of the state and chemistry of the IGM in the first billions years of the Universe 
- Ahn et al. 2015; Subrahmanyan et al. 2015 

• Unique studies (i.e. CD) beyond the Epoch of Reionization which will remain out of reach 
of most of not all other (planned) facilities in particular of the first (Pop-Ill stars and X-ray 
heating sources) - Mesinger et al. 2015; Semelin et al. 2015 

• Probing the impact of bulk-flows during the later Dark Ages and the Cosmic Dawn allowing 
physics during the Dark Ages and CMB to be probed in a unique manner - Maio et al. 2015 

• Strong synergy with intensity mapping of the CO, CII, Lyman-a lines, possible molecular 
lines from primordial collapsing gas (e.g. H 2 and HD), as well as the kS-Z effect and NIR/X- 
ray mission, as well as with many other (planned) facilities in space and on the ground — 
Chang et al. 2015; Jelic et al. 2015 

We conclude to say that direct detection of neutral hydrogen at high redshift, beyond the statistical 
approaches of current experiments, as well as probing the Cosmic Dawn will be defining abilities 
of the SKA (1 and 2). SKA can probe the universe back to 100 Myr after the Big Bang, well beyond 
even the limitations of JWST 9 . This unique science can be carried out in a three-tiered survey of 
7500 hrs or total on-sky time. 
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